Clarifying spherical collapse in coupled dark energy cosmologies 
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The spherical collapse model is often used to follow the evolution of overdensities into the non- 
linear regime. We describe the correct approach to be used in coupled dark energy cosmologies, 
where a fifth force, different from gravity and mediated by the dark energy scalar field, influences 
the collapse. We reformulate the spherical collapse description by deriving it directly from the set 
of nonlinear hydrodynamical Navier-Stokes equations. By comparing with the corresponding rel- 
ativistic equations, we show how the fifth force should be taken into account within the spherical 
collapse picture and clarify the problems arising when an inhomogeneous scalar field is considered 
within a spherical collapse picture. We then apply our method to the case of coupled quintessence, 
where the fifth force acts among cold dark matter particles, and to growing neutrino quintessence, 
where the fifth force acts between neutrinos. Furthermore, we review this method within standard 
cosmologies and apply our analysis to minimally coupled quintessence. We also check past results 
for early dark energy parametrizations. 
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I. INTRODUCTION 

A wide variety of theoretical cosmological models can 
be challenged and discriminated thanks to predictions 
on structure formation. At the nonlinear level the be- 
havior of ACDM cosmologies, in vifhich the role of dark 
energy (DE) is played by a cosmological constant, can 
significantly differ from dynamical dark energy models. 
In more realistic scenarios, these allow for DE couplings 
to other species. Interacting dark energy cosmologies in- 
clude: coupled quintessence (DE evolution is coupled to 
dark matter) |l MlOl|: g rowing neutrino cosmologies |ll| - 
[l3| and MaVaNs ([lj4iil and references therein) (DE 
is interacting with neutrinos); so-called modified gravity 
theories such as scalar-tensor theories, including F(R) 
and extended quintessence |2CH29f . In all these cosmolo- 
gies a fifth force is present, acting on species whose evo- 
lution is coupled to the DE evolution. The presence of a 
fifth force, mediated by the DE scalar field (the cosmon, 
seen as the mediator of a cosmological interaction) can 
modify structure formation in a significant way |30l . |3l| , 
in particular at large scales [32| . In view of future data, it 
is therefore important to understand how these theories 
behave when density perturbations reach nonlinearity. 

While up to now N-body simulations represent the 
best way to numerically evolve structures, other semi- 
analytical methods have been used to follow pertur- 
bations into the nonlinear regime, either using spher- 
ical collapse |33l438| or other alternative methods [39|- 
4l|. In particular, spherical collapse has been used in 
several occasions in literature for ACDM 1351 l36l |42| . 
minimally coupled quintessence models [43145 1| . coupled 
quintessence 52, |53| and when parametrizing early dark 
energy contributions 54l456l| . 

In this paper we give a detailed description of the 
spherical collapse method and clarify some tricky issues 



in its applications. We lay particular focus on the cal- 
culation of the extrapolated linear density contrast at 
collapse 5c, a quantity of major interest within a spher- 
ical collapse description, often used in a Prcss-Schechter 
[57| approach to estimate dark matter halo mass distri- 
butions. 

After reviewing results for standard cosmologies like 
ACDM, we consider the case in which a fifth force is 
present in addition to standard gravitational attraction, 
as in the case of all the interacting dark energy models 
mentioned above. The inclusion of the fifth force within 
the spherical collapse picture requires particular atten- 
tion. Spherical collapse is intrinsically based on gravita- 
tional attraction only and cannot account for other exter- 
nal forces unless it is suitably modified. The dynamics in 
the spherical collapse model are governed by Friedmann 
equations. Hence, only gravitational forces determine the 
evolution of the different scale factors and, in turn, of the 
density contrast. 

A detailed comparison between the linearized spheri- 
cal collapse picture and the linear relativistic equations 
allows us to first identify the presence (or absence), in 
the spherical collapse picture, of terms which are a direct 
signature of the coupling already at the linear level. We 
show how spherical collapse necessitates to be suitably 
modified whenever an additional force other than gravity 
is present and is big enough to influence the collapse. We 
use this comparison also to show that a standard treat- 
ment of spherical collapse may lead to problems even in 
the uncoupled case when treating inhomogcneities in the 
scalar field. 

A modification of the spherical collapse picture is in- 
deed possible, via a nonlinear analysis of the model. We 
derive the set of second order differential equations for 
the density contrast from the nonlinear Navier-Stokes 
equations described in j32l | , extending an idea from [58| . 



We show how Sc can be evaluated directly from these 
equations and how they can serve as a starting point for 
a reformulation of spherical collapse. Our results match 
the numerical solution of the nonlinear hydrodynamical 
equations solved as described in |32|.|58|. 

We apply our method to coupled quintessence scenar- 
ios where a coupling is present among dark matter par- 
ticles, comparing our results with alternative methods 
presented in the past [52l [53[ . As a further application 
of our method, we consider for the first time spherical 
collapse within growing neutrino models, where an inter- 
action is active among neutrinos: in this case we obtain 
an extrapolated linear density at collapse Sc which shows 
an oscillating behavior, a characteristic feature of the in- 
teraction. 

Finally we confirm results found in [56| on spherical 
collapse and early dark energy (EDE). 

In Sec. [n] we recall the spherical collapse model and 
its applications to standard cosmologies (Sec. Ill A[) . In 
Sec. mil we focus on spherical collapse in presence of a 
fifth force, taking the case of coupled quintessence as an 
example of fifth force cosmologies. In Sec. lIVI we demon- 
strate that the standard spherical collapse leads to wrong 
results when applied to coupled quintessence: indeed, 
by comparing with full rclativistic equations (Scc lIV A]) 
we demontratc that the spherical collapse equations lack 
terms that arc essential in the presence of a fifth force 
fSec lIVBJ) and can lead to incorrect results when an inho- 
mogeneous DE scalar field is included within this frame- 
work (Sec. lIVCp . Consequently, in Sec.|Vl we illustrate 
how the spherical collapse can be correctly reformulated 
in coupled scenarios by basing it on the full nonlinear 
Navicr-Stokcs equations for the respective model. We 
further comment on the careful choice of initial condi- 
tions in Sec. lVAl We apply the derived formalism to give 
results for coupled quintessence (Sec. IV Bp and growing 
neutrinos (Scc. lVIj) . Finally, we use the described frame- 
work to confirm results found in [56J for uncoupled early 
dark energy. 



II. SPHERICAL COLLAPSE 

Consider a cold dark matter density perturbation 
within a homogeneous background Universe. Under the 
effect of gravitational attraction the perturbation grows, 
possibly entering the nonlinear regime, depending on the 
scale of the perturbation. A popular method often used 
to follow the evolution of cold dark matter (CDM) struc- 
tures during the first stages of the nonlinear regime is 
the spherical collapse model. In its original applications 
J33l436l |. it is assumed that the initial overdensity obeys 
a top hat profile 



where r{t) specifies the radius of the top hat and s is 
the spherical coordinate indicating the distance from the 
center of the perturbation. 8(r(i) — s) is the top hat 
function, equal to 1 for s < r{t) and otherwise. The 
amplitude of the top hat is given by Spo and is evolving in 
time. As a consequence of Birkhoff 's theorem of General 
Relativity, which ensures that the dynamics of the ra- 
dius r{t) are governed only by the enclosed mass, the top 
hat "bubble" is conveniently described as a closed Uni- 
verse where the total density p = pcrit + Spm exceeds the 
critical density pcrit due to the presence of CDM density 
perturbation. 

Hence, all densities and geometric quantities are 
treated according to the Friedmann equations: 



H' 



1 v^ K 

a 



(2) 
(3) 



Here the "scale factor" is given by the radius of the bub- 
ble r{t), commonly normalized to match the background 
scale factor a{ti) at some initial time ti. The correspond- 
ing Hubble function of the bubble is indicated by H . 
Eq. ([2]) explicitly contains a curvature term K\ Eq. ([3]), 
albeit the lack of an explicit curvature term, is still de- 
scribing a closed Universe as the sum of the densities 
on the right hand side exceeds the critical one. Note 
that throughout this work densities have been normal- 
ized in units of the square of the reduced Planck mass 
Af2 = (87rGAr)"\ 

The bubble is embedded in a homogeneous Fricdmann- 
Robertson- Walker (FRW) background characterized by 
a scale factor a{t) and a corresponding Hubble function 
H = a/a. We use a bar to indicate background quanti- 
ties. For clarity, we recall the Friedmann equations de- 
scribing the homogeneous and flat background Universe: 



H 



= 3 H ^" ' 

Q. 



a 6 



(4) 
(5) 



5p,^{t,s) = 6po{t)Q{r{t)-s) 



(1) 



Note that throughout this work we neglect baryonic 
components. For simpler notation we refer to CDM by a 
subscript m. 



A. Applications to standard cosmologies 

Spherical collapse can be safely applied to the case of 
Einstein de Sitter (EdS) cosmologies (in which r2„i = 1), 
and to ACDM models. In this case, the energy density 
of matter pm, appearing on the right hand side of Eq. ^ 



and ^, is conserved both inside and outside the over- 
density: 



Pm = -SH{l+Wm)Pn 
Pm = -3H(1 + Wm)pn 



(6) 
(7) 



The nonhnear density contrast is defined by 1 + 5m = 
Pm/ Pm and is determined by the above equations. The 
hnear density contrast evolves according to well known 
linear perturbation theory J59l |60| and satisfies the linear 
equation: 



S. 



'rn,L 



i 3 , 

'^HSmX — -^H''Vtra5m,L 







(8) 



Equations ([3]) - (H)) can be integrated numerically. We 
start the integration at some initial time tin in which the 
total energy density in the bubble is higher than the crit- 
ical energy density, due to the presence of the CDM over- 
density 5m- Equation ([3|) provides r(z), which is shown 
in FigHlfor a ACDM model with I^a = 0.7 and for three 
different overdensities (Ji„ = 1 • 10"^, 2 • 10"^, 3 • 10"^ 
for Zin = 10''): r{z) first increases as the bubble expands 
with the background; then, it reaches a maximum value 
(turnaround) in which comoving velocities become zero; 
finally, the bubble collapses, the radius tends to zero and 
the nonlinear density contrast 5m increases rapidly. The 
redshift of collapse depends on the amplitude of the ini- 
tial perturbation. The higher this is, the earlier the over- 
dense region will collapse. The corresponding value of 
the linear density contrast extrapolated at the time of 
collapse is usually referred to as 5c and represents one of 
the key ingredients for a Press-Schechter analysis, which 
gives statistical estimates of the cluster distribution in 
space. We will not go into the Press-Schechter procedure 
here; instead, we will focus on the calculation oi5c within 
the spherical collapse analysis, putting in evidence how 
this calculation has to be carefully performed depending 
on the underlying theoretical model. 

In an Einstein de Sitter model the linear density con- 
trast at collapse can be calculated analytically [35|, |36[ : 
it is equal to a constant value independent of the redshift 
of collapse Zc 



5c = (3/20) (127r) 



2/3 



1.686 



(9) 



Note that we define Zc as the redshift at which r — > 0. In 
a ACDM model one expects this value to decrease for late 
collapse times |35l . |36| , when dark energy dominates over 
matter and leads to cosmic acceleration, slowing down 
structure formation. In Fig. [5] we plot 5c{zc) for both 
EdS and ACDM. 

It is also common to analyze the (nonlinear) density 
contrast at virialization. From the virial theorem one 
may deduce [35|, ISfl |42|, |46| that a given bubble virializes 
whenever it has collapsed to half its turnaround radius. 
In an EdS Universe, the density contrast at virialization 
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FIG. 1. Evolution of radial parameter r{z) for different initial 
overdensities in a ACDM model with Q.m = 0.3, Q.K = 0.7, as 
used throughout this work. For comparison, we have included 
the background scale factor a (dotted, pink). 
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FIG. 2. Extrapolated linear density contrast at collapse 5c vs. 
redshift at collapse Zc for a ACDM (solid, red) and an EdS 
(double-dashed, black) model. 



For the ACDM model, an increase is observed for late 
collapse times. This corresponds to the fact that, in pres- 
ence of dark energy, it takes longer for structures to viri- 
alize, with a corresponding higher value of (Svir, as shown 
in FigH 

To better illustrate the effect of A on both (5vir and i5c, 
we have plotted 5m,L and 5m, nl in FiglH as well as the 
radius r for a fixed initial overdensity in both an EdS 
and a ACDM Universe. It can be seen that the later 
virialization in ACDM leads to an increase in 5vir- On 
the other hand, the smaller linear growth rate reduces 
the extrapolated linear density contrast 5c in ACDM. 
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[63| . Viable models in this direction often involve the 
presence of a coupling between the dark energy scalar 
field, referred to as "cosmon" or "quintessence", and 
other components in the Universe such as cold dark mat- 
ter [l|, \2^ or neutrinos lll-ll4 |32| . The presence of an 
interaction that couples the cosmon dynamics to another 
species introduces a new force. This "fifth force" is act- 
ing between particles (CDM or neutrinos in the examples 
mentioned) and is mediated by dark energy fluctuations. 
Whenever such a coupling is existent, spherical collapse, 
whose concept is based on gravitational attraction, has 
to be suitably modified. In the following sections we will 
present some examples of quintessence models in pres- 
ence of a fifth force and show how the latter can be taken 
into account. 



FIG. 3. Nonlinear density contrast at virialization for a 
ACDM (solid, red) and an EdS (double-dashed, black) model. 



III. COUPLED QUINTESSENCE 
COSMOLOGIES 
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FIG. 4. Linear and nonlinear density contrasts in ACDM 
(red) and EdS (green) models, as well as the corresponding 
radius functions (blue, short-dashed and pink, dotted, respec- 
tively). The upper double-dashed black line marks the turn 
around radius rta, the lower one rta/2. The later virialization 
time in ACDM (rta/2 is reached substantially later, when the 
blue dotted line and the lower black double-dashed lines in- 
tersect) leads to an increase of S^iv The overdensities collapse 
when the radii go to zero and 5c is given by the value reached 
by the linear curves at this redshift. Although collapse hap- 
pens later for ACDM, the linear growth is suppressed at late 
times (red dashed line in comparison to green dashed line). 
This leads to a decrease in 5c ■ 



Alternatively to a cosmological constant, dark energy 
can be described by a dynamical energy component, such 
as a. quintessence scalar field rolling down a potential 



6ll |62| . A meaningful quintessence model should natu- 



rally explain why dark energy dominates over cold dark 
matter only at recent times; this happens to be difficult 
to achieve within minimally coupled quintessence mod- 
els, which are often fine-tuned as much as a ACDM model 



The first set of cosmologies in presence of a fifth force 
that we consider is coupled quintessence (CQ): here the 
evolution of the quintessence scalar field (from hereon we 
refer to it as the 



cosmon 
The cosmon (t> 



611) is coupled to CDM [i 



|6|, |22|, |30[ • The cosmon interacts with CDM particles 
whose mass m{(j)) changes with 0. This set of cosmologies 
is described by the Lagrangian: 

C = -^d^4>d^4> - Ui^) - m{4>)^P^ + ^kinM , (10) 

in which the mass of matter fields -0 coupled to DE is a 
function of the scalar field cj). 

The homogeneous fiat background follows the set of 
equations described in [l|, [ij] . The Universe evolves in 
time according to the Friedmann and acceleration equa- 
tions: 



H' = 



gE' 



and 



Q^[pM+5Wa)] 



(11) 



(12) 



where the sum is taken over all components a in the Uni- 
verse. A crucial ingredient is the dependence of CDM 
mass on the cosmon field (f), as encoded in the dimen- 
sionless cosmon-CDM coupling /3, 



/? = - 



din 7 



(13) 



For increasing (f) and 
decreases with time 



> the mass of CDM particles 



m 



me 



-P4, 



(14) 



where m is a constant and /? is also fixed to be a con- 
stant in the simplest coupling case. The cosmon field 
(j) is normalized in units of the reduced Planck mass 
M = (87rGiv)~^^^, and P ^ 1 corresponds to a cosmon- 
mediated interaction for CDM particles of roughly grav- 
itational strength. 

For a given cosmological model with a set time depen- 
dence of (p, one can determine the time evolution of the 
mass m(i). The dynamics of the cosmon can be inferred 
from the Klein Gordon equation, now including an extra 
source due to the coupling to CDM: 



3h} +^=l3p„ 



(15) 



We choose an exponential potential |6ll.l62l.l64l.l65|: 



V{(j)) = M^U{<I>) ^ M^e- 



(16) 



where the constant a is one of the free parameters of our 
model. Note that our analysis, however, is more gen- 
eral and can be applied in presence of any quintessence 
potential. 

The homogeneous energy density and pressure of the 
scalar field (f) are defined in the usual way as 



-C/(0) , p0 



^U{(f)) , w^ 



p^ 



(17) 



Finally, we can express the conservation equations for 
dark energy and coupled matter as follows [ij, 122] : 



P<t> 



-3H{1 



W^)P4, + P4>Pn 



(18) 



The sum of the energy momentum tensors for CDM and 
the cosmon is conserved, but not the separate parts. We 
neglect a possible cosmon coupling to baryons (b) or neu- 
trinos {v), so that pb^i, = —3H{1 + 'Wb.i^)pbM- 

For a given potential p6|) the evolution equations for 
the different species can be numerically integrated, giving 
the background evolution shown in FiglS] (for constant 
/3). For a detailed description of attractor solutions in 
this context see [l|, \2^ |j 



It will later prove useful to understand the evolution of 
perturbations within coupled models in the linear regime. 
The relativistic calculation in coupled quintessence is de- 
scribed in detail in [J, llSl, |32|. Here we just recall the re- 
sulting second order equations (in Fourier space) for Sm,L 
and the perturbation of the scalar field 5(f> in Newtonian 
gauge (in which the nondiagonal metric perturbations are 
fixed to zero): 



'Jm,L 



-2H iSrn,L 

-13 CJ) {5m,L 

-3<P-(3S4> 



-/3i 



+ (fc/a)2 ($ + /?, 



3$(2i/ 



P$) 



(19) 
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FIG. 5. Energy densities of cold dark matter (solid), dark 
energy (dotted) and radiation (long dashed) are plotted vs 
redshift. We take a constant jS — 0.1, with a = 0.1. 



= -3HSq) - t/,^^ S(f> + I5p^ ((S„,,L + 2$) 

-(fc/a)2j0-2$C/,^^ + 40$ . (20) 



Since the spherical collapse is intended to model the 
nonlinear evolution in the Newtonian limit, we are inter- 
ested in the case in which k ^ aH. Following [4|, |6|, we 
obtain 

^™ = (/30-25)5™ + (fc/a)2(l + 2/32)$ , (21) 



eScj^^Pa'-p^S^ . (22) 

and the gravitational potential is approximately given by 

fc2$ ^ 1^2 ^ -^^^ ^ (23) 

where we have assumed that no anisotropic stress is 
present, so that $ = — ^I*. We can then define an effective 
gravitational potential as 



$eff = $ + PH 



(24) 



In real space (comoving spatial coordinates) and after 
substituting the expressions for $ [Eq. ([25]) ] and for (50 
[Eq. (P^ ]. we get the modified Poisson equation: 



A$, 



-Pm^m (1 + 2/?2) 



(25) 



Cold dark matter then feels an effective gravitational con- 
stant 



Gcff = GA.[l + 2/32] , 
where Gn is the usual Newton's constant. 



(26) 



The first term on the right hand side of Ea. (pT|) in- 
cludes the expansion damping, modified by the velocity 
dependent term /3(/), which accounts for momentum con- 
servation; the last term on the right hand side specifics 
the presence of the fifth force. 



IV. STANDARD SPHERICAL COLLAPSE AND 
COUPLED QUINTESSENCE 

We will now apply the framework described in Sec. 
im as it is to CQ. This approach has been used, for ex- 
ample, in j52l |. We will show, by comparison with the 
perturbation equations recalled in SeclIIH that this does 
not correctly model the evolution of nonlinear structures 
in coupled quintessence. 

For this purpose, consider the standard spherical col- 
lapse equations ([2]) and ([3]), in which the densities on the 
right hand side satisfy the coupled conservation equa- 
tions: 



Pr = -AHpr + Tr , 
Pm = -iHpm - Pi'Pm , 
/90 = -3H{p^ + p^) + fi(j)Pcdm + r^ 

or equivalently the Klein-Gordon equation: 



(27) 
(28) 
(29) 



(30) 



Here additional source terms F,, and F^ may account for 
possible differences between the bubble and the back- 
ground coraponents for radiation and the scalar field, re- 
spectively [4514471 . [52 ] . In case of clustering dark energy 
and/or radiation, both source terms are set to zero: 



F^ =0 
F^ = 



(31) 
(32) 



If both radiation and the scalar field are to be homoge- 
nous, i.e. they behave in the bubble as in the background, 
the source terms are defined as: 

Tr=A{H~H)p,. , (33) 

T^ = i{H-H){p^~p^)+(i{^^$) . (34) 

In order to account for a fractional outflow of dark energy 
or radiation, one may suitably interpolate between the 
two values |45l - l47l .[5 ^ . 



A. Comparison with relativistic equations 

We will now show that the approach described by Eqs. 
© and dH) together with ([27] - [M]) is incorrect. The 
actual fifth force term is entirely missing from the equa- 
tions. 



As a starting point, we remark that Eq. p8|) for the 
background CDM density and dMl) for the bubble CDM 
density can be directly integrated [l| to yield 



Pn 



J4>i. 



./9 0i„ 



-3 



.-/3 



-P<t> 



(35) 
(36) 



The density contrast 5m is then given by 
1 



r Pm 

Orn = ^— 



[l + Sra^rde''"- — 
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- 1 e 

r 



-/3 50 



where we have introduced S(j) = 4 
second time derivatives of S„i read 



(37) 
The first and 



Srn^3{l + S.m){H-H)-l3S4>{l + Sm) , (38) 

1 + dm ^ ' 

(39) 



\ + bm ^ " V" "I 

-/3<5^(l-f,5™) . 

We can substitute {U — H) and (H — Hj using Eqs. ^ 

and ([3]). Taking the square of Eq. ([55| and inserting it 
into Eq. ([5^ we obtain: 



-2Hi^6m+/3S^il+Sm) 
- (1 + Sm) ^ {Spa + 36pa) 



3 1 



> 



> 



^(364>il + Sm) 



(40) 



This is the evolution equation for the density contrast Sm, 
as derived directly from spherical collapse applied to cou- 
pled quintessence. Usually, one considers only cold com- 
ponents to actually cluster, reducing the sum in Eq. ()40|) 
to one over CDM only. For the moment we still allow 
for an inhomogenous scalar field and therefore also for 
nonvanishing Sp^ and Sp^. 

The evolution of S(f) = 4>—(f> is determined by combining 
the Klein-Gordon equations for the bubble ((30|) with that 
of the FRW background P^ : 



(41) 




- c/.^L- [/., 



l3Sp„ 



Again, we remark that Eqs. (pH)) and (PT|) are obtained 
by applying standard spherical collapse equations (jH]) and 
([5|) to coupled quintessence simply by adding a coupling 
in the conservation equations ((?7)) - (|34p . If linearized. 



Eqs. (|40)) and (|4ip read as shown in the left column of 
Table H 

In Table H] we compare the equations found for Sm and 
S(j> obtained from standard spherical collapse to Eqs. p^ 
and (POI obtained from the fully relativistic theory in 
Newtonian gauge and shown in the central column of 
Table |T1 In the right column we display the relativistic 
equations within the Newtonian limit, corresponding to 
Eqs. (|22[l and pTj) . Some remarkable problems become 
evident. 

Note that we have chosen to display the relatistic equa- 
tions in Newtonian gauge to analyze time derivatives of 
64> in the spherical collapse equations, as well as to weigh 
the importance of different terms in the Klein-Gordon 
equation when going to small scales. One should keep 
in mind that at large scales, these equations are gauge 
dependent. 



B. Lack of the fifth force in spherical collapse 

As compared to both the Newtonian and relativistic 
equations, major terms are missing in the standard spher- 
ical collapse scenario: 

• Comparing (a) to (b) and (c) in Table U no term 
proportional to /?(/> appears. Depending on the 
strength of the coupling f3, this term, originating in 
momentum conservation, can be of great relevance. 
For /3 ~ 1, it can significantly alter structure for- 
mation when correctly considered in the vectorial 
velocity equations, as shown in [30|. For large cou- 
plings, as e.g. in a growing neutrino scenario we 
will discuss later, it is less important, since the cos- 
mon (j) is almost constant at late times. Comparing 
(g) to (h) and (i), one notices a sign reversal in 
front of a friction-like term. This will also lead to 
wrong results. 

• Comparing (a) to (b) and (c), as well as (g) to (h), 
in Table U terms proportional to fik^Scj) are absent 
in the spherical collapse equations; this term is ex- 
actly what provides the fifth force. Its absence in 
(g) leads to a sign reversal of /3^ Spm,L as compared 
to (i), thus yielding an incorrect effective Newton's 
constant. 

• Comparing (d) to (e) and (f) in Table U terms 
proportional to k'^Scj) are absent in the spherical 
collapse equations; remarkably, this term does not 
depend on P and is therefore missing even in the 
uncoupled quintessence scenario. 

The lack of terms proportional to (3k'^S4> leads to a de- 
scription which does not correspond to the desired cou- 
pled quintessence scenario: indeed, these are exactly 
the terms responsible for the fifth force, originating in 



m 

Eq. 



and leading to an efi^ective gravitational force as in 
26|) . In other words, we point out that the standard 

52| does not 



spherical collapse, as used for example in 
include the main ingredient of coupled quintessence. A 
fifth attractive force acting between CDM particles and 
mediated by the cosmon is absent, although densities are 
indeed coupled to each other as in pS]) - (pO)) . The rea- 
son for this can be seen as follows: spherical collapse 
is by construction based on gravitational dynamics and 
cannot account for other external forces unless appropri- 
ately modified. The dynamics in the spherical collapse 
models are governed by the usual Eriedmann equations, 
which arc particular formulations of Einstein's field equa- 
tions. Hence, only gravitational forces determine the evo- 
lution of the different scale factors and, in turn, of the 
density contrast. We note that, though in the limit of 
small couplings the difference can be small, for strongly 
coupled scenarios a completely different evolution is ob- 
tained. This is simply connected to the fact that for small 
couplings gravity is still the crucial ingredient to fuel the 
collapse. 



C. Inhomogeneity of the scalar field 

The issue of whether the scalar field should be consid- 
ered to be homogeneous (with the cosmon inside the top 
hat given by the homogeneous background field) or not 
has also been addressed in literature. In particular, one 
could try to compare a homogenous scalar field (j) to an 
inhomogeneous one by appropriately fixing F^ in Eq.()30|) 
to the expression p4p an d ([5^ respectively. This com- 
parison led, for example, [53 to find differences between 
the homogenous and inhomogeneous cases. 

The difference found following such procedure is, how- 
ever, not caused by the fifth force, which, as shown, is 
not present. Even without any fifth force or coupling, we 
have noticed in Table HI comparing (d) to (e) , that terms 
proportional to k'^6(f> are absent in the spherical collapse 
equations. 

Evaluating the clustering Scj) using a spherical collapse 
scenario as given by (d) in Table U leads to effects which 
do not correspond to the relativistic behavior. In fact, 
in absence of the term —k'^Scj), spherical collapse overes- 
timates the time dependence of the scalar field perturba- 
tions as soon as S(j) is assumed to be different from zero: 
if, for example, 54>in > initially, the S(f> obtained from 
the bottom left equation (d) in Table U is bigger than it 
would be if the term —k'^S(j) appearing in the relativistic 
equations was actually present. Hence, within spherical 
collapse, all time derivatives of the cosmon are overesti- 
mated and, as a consequence, Sm is incorrectly reduced. 

We remark that this reasoning also applies to or- 
dinary, uncoupled quintessence. Here the question of 
(in)homogeneities in the scalar field was addressed in var- 
ious works, e.g. [45|-|49|. Also in this case, one may gain 
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Relativistic perturbations 
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+/3^4</' + t^,00/35</' 
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(g) 


(h) 


(i) 





TABLE I. Comparison between linearized spherical collapse and fully relativistic linear evolution equations. The third row is 
a combination of the first two rows. 



some insight by considering the equations in Table U now 
for /3 = 0. In the relativistic description, scalar field per- 
turbations will decay due to the presence of the term 

— (k/a)'^5<j>, until the latter may eventually be countered 
by gravitational contributions. In the spherical collapse, 
however, this decay is only driven by the scalar mass 
term —U^4,4,S(j). For a light scalar field and/or sufficiently 
small scales, this term is smaller than the missing term 

— {k/a)'^S4> and therefore scalar field inhomogeneities are 
incorrectly overestimated. The error may be substan- 
tially reduced if heavy scalar fields are considered. In 
this case, the presence of the large mass term V"{(p)5(j) 
in the perturbed Klein-Gordon equation may make up 
for the lack of spatial gradients. 

In conclusion, we have shown that applying the spher- 
ical collapse equations to coupled quintessence by merely 
modifying the conservation equations can lead to results 
which do not correspond to the wanted cosmological sce- 



nario; this procedure in fact not describe the nonlinear 
evolution in CQ. It is, however, possible to amend the 
above model to properly include the fifth force whenever 
a coupling is present. In order to illustrate and justify 
that, we consider the nonlinear hydrodynamical evolu- 
tion equations within coupled quintessence scenarios. 



V. HYDRODYNAMICAL SPHERICAL 

COLLAPSE: A CONSISTENT APPROACH TO 

COUPLED QUINTESSENCE 

In uncoupled, purely gravitational cosmologies the 
spherical collapse can be derived from the hydrodynami- 
cal Navier-Stokes equations. This is a consequence of the 
fact that the Friedmann equations in presence of non- 
relativistic components can be derived from Newtonian 
gravity. We will demonstrate now that this is also possi- 
ble in the presence of external forces, basing our analysis 



on an idea first developed in [58|. 

In order to derive the correct formulation in coupled 
quintessence, we consider the full nonlinear evolution 
equations in coupled cosmologies within the Newtonian 
limit: 



(5m = 


-Vm V(5„i - (1 + (5m) V • Vm 


(42) 


Vm = 


-(2if-/3 0)v„, -(v,„V)v™ 






-a~2V($-/3(5(?i) 


(43) 


AS(j) = 


"/? a^ (5pm 


(44) 


A$ = 


— :r ^ '(5/9„ 


(45) 



These equations can be derived both from the nonrel- 
ativistic Navier-Stokes equations and from the Bianchi 
identities in the appropriate limit in presence of an ex- 
ternal source 1591 . 



v,t; 



Q^ = -md^cj, 



(46) 



where TJ is the stress energy tensor of the dark matter 
fluid. They are valid for arbitrary quintessence potentials 
as long as the scalar field is sufficiently light, i.e. mi^Scj) = 
V"{(f>)S(j> <C A5(f> for the scales under consideration. For 



a more detailed discussion of the equations, see [32 



We are working in comoving spatial coordinates x and 
cosmic time t. The sign in Eq. ([15|) was chosen to match 
Eq. ((25)) . Note that v^ is the comoving velocity, related 



to the peculiar velocities by v^ = Vpec/a- The sum in 
Eq. (|45|) is to be taken over all clustering components; as 
an important consequence of the Newtonian limit, the 
cosmon is explicitly excluded. 

In order to obtain a correct description of the spherical 
collapse model, we are interested in the evolution of a 
top hat, spherically symmetric around x = 0. We note 
that the below derivation is not limited to a top hat but 
holds for the amplitude at x = for generic spherically 
symmetric profiles. From simple symmetric arguments 
we may infer 



V(5„ 



.(0,t) = 



which changes 



to 



x=0 



[(1 + (5„) V • V, 



mj Tr=i 



x=0 



(47) 



(48) 



We now want to relate Eqs. (|42p - ([l5)) to the spherical 
infall: it is therefore useful to combine ((42]) and (|43| to 
give a second order equation for S„n taken at x = 



[(l + ,5„)V-v„]Uo , (49) 



c=0 1 + 5n 



where wc have used (|47p and (|^5)) . Inserting the diver- 



gence of (|43l) yields 



x=0 



-i2H-/3^)6„ 



1 



A$, 



1 + (5„, ( 

+ (l + (5„,)V(v„V)v„|^^o 



(50) 



Note that ^^s is defined as in (|24| and obeys the Laplace 
equation (PS|) , as can be seen by combining (^^ and (|45p . 
The first three terms in ([50]) can be evaluated straight- 
forwardly at X = 0. To rewrite the last term we use the 
identity 



V(vm V)v,i 



lx=0 



(V ■ v„ 
SI 



lx=0 



(51) 



3 (l + (5,„)2 



which holds for spherically symmetric situations and is 
rederived in the Appendix. Inserting this into expression 
([76|) and subsequently into ((50)) yields the final expres- 
sion for the evolution of the top hat density amplitude 
(writing 6 instead of (5|x=o) 




Linearization leads to: 




which corresponds to the relativistic equation ((21)). Here 
we recall that the effective gravitational potential, given 
by ((24)) , follows the modified Poisson equation p5p which 
we rewrite here for convenience: 



A$eff=-yPm^m(l + 2/32) 



(54) 



Equations (|52p and (|53p are the two main equations 
which correctly describe the nonlinear and linear evo- 
lution for a coupled dark energy model. They describe 
the dynamics of a spherical top hat as it follows from 
relativistic perturbation theory in the Newtonian regime 
and they can be used, among other things, for estimating 
the extrapolated linear density contrast at collapse 6c in 
the presence of a fifth force. To our knowledge it is the 
first time that the second order equations ((5^ and ((55)) 
are presented in this way. 
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We will now demonstrate that we may easily reformu- 
late Eqs. (|52p and (|53p into an effective spherical col- 
lapse: we can combine tlicni to derive an equation for 
the radius r which extends Eq.Q to the case of coupled 
dark energy. To do so, we consider a spherical bubble of 
radius r containing the CDM overdensity dm- Particle 
number conservation yields 



(1 



^n, 771. my 



r / 



(55) 



where n is the number density of CDM particles and 
5n = Sn/n. 

We demand the scale factors r and a to be equal ini- 
tially, i.e. Oin = Tin. Further, we assume that the mass 
of CDM particles is the same inside the bubble and in 
the background. Note that this is not a limitation, but 
merely a prescription that we have employed in order to 
obtain an equation for the scale factor r in a form which 
is analogous to the original Friedmann equation ^ . We 
obtain 



1 + 5,n = (1 + (5,n,i„) - 

\r 



(56) 



The first and second time derivatives of (5„, then read 



<j,„ = 3 (1 + 5„0 I - - - 
\ a r 



<5m = 3(1 + S,n) 



(57) 




l + Sra ' 

which we can combine appropriately to yield 
- • 4 (5^ / n r 

5m = ^2H5m + ^-^^+Kl + 5m) 

3 1 + dm \ar 



(58) 



(59) 



Comparison to ([5^ and insertion of the background 
Friedmann equation ([2]) gives the evolution equation for 
the bubble radius 



f 
r 


\ 

-\p^5pm 


-0- 




(60) 



Equation (pO)) . equivalent to the one used in [53|, de- 
scribes the general evolution of the radius of a spherical 
overdense region within coupled quintessence. Compar- 
ing with the Friedmann equation ([3]) we notice the pres- 
ence of two additional terms: a friction term and the cou- 
pling term /3^ 5pm', the latter is precisely the term respon- 
sible for the additional attractive fifth force. Note that 



the "friction" term is velocity dependent and its effects 
on collapse depend, more realistically, on the direction 
of the velocity [30[, information which is not contained 
within a spherical collapse picture. 

We conclude that one may indeed apply the spherical 
collapse model to coupled dark energy scenarios. How- 
ever, it is crucial to include the additional force term in 
the equations. 

Note that the outlined procedure can easily be gen- 
eralized to include uncoupled components, for example 
baryons. In this case, the corresponding evolution equa- 
tion for (5fc; will be fed by $off = ^- This yields an evo- 
lution equation for the uncoupled scale factor 7'„c that is 
equivalent to the regular Friedmann equation ([3]). 



A. Methods and initial conditions 

To provide maximum stability and to rule out a de- 
pendence on initial conditions, we directly integrate Eqs. 
(|52p and (|53p for the nonlinear and linear density con- 
trasts, together with the corresponding background equa- 
tions and the Klein-Gordon equation ([T5]) for the scalar 
field. The radial parameter r[z) may equivalently be ob- 
tained by integrating Eq. ([60)) or by directly applying the 
relation (|56p. The following initial conditions at the ini- 
tial redshift Zjn were imposed: 

* ^m,in — ^m,L, in 

• 5m.L.\n = 3(1 -I- 5m,LM){H\n " ^in) = 0, as initially 
the Hubble functions of background and overden- 
sity evaluate to the same value. 

The value of the extrapolated linear density contrast at 
collapse 5c can be obtained by stopping the evolution of 
Eg. ([53)) when 5m as obtained from ([52)) goes to infinity, 
i.e. the overdensity collapses. If we then vary the initial 
conditions, leading to different collapse redshifts Zc, we 
arrive at a redshift dependent expression for this critical 
density, 5c = 5c{zc)- Equivalently, one may vary the 
initial redshift Zjn, keeping 5m:vi fixed. To be sure of 
starting the integration when densities are still linear, 
we find that it is necessary to work in a range of initial 
ovcrdcnsities with (5m. in < 10^'^. 



B. Results 

We depict the evolution of (5m(z) and (5m,L(^) for differ- 
ent initial redshifts in Fig. [Hj For this plot, we used sam- 
ple parameters a = 0.1 and /3 = 0.1. We have also plot- 
ted the linear density contrast at collapse 5c{zc) for three 
coupled quintessence models with a = 0.1 and p = 0.05, 
0.1 and 0.15 in Fig. [7| We note that these as well as 
all subsequent results are valid under the hypothesis in 
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which the linear extrapolation traces the nonhnear be- 
havior when a fifth force is present. 




1000 



1+z 



FIG. 6. CDM linear and nonlinear perturbations for different 
initial conditions. 



1.71 




for this increase is quite simple. In Eas. ((52|) and ((53)) two 
terms lead to an enhanced growth: the fifth force term in 
the effective potential and the reduction of the damping 
~P4>. In the linear equation, they are always of compa- 
rable strength. In Eg. ([5^ . however, the damping will be 
negligible once (5m '^ 1 as it only enters the equation lin- 
early. The enhancement of growth is then weaker than 
in the linear equation and 5c grows with increasing /3. 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 
P 



FIG. 8. Extrapolated linear density contrast at collapse 
5c for coupled quintessence models as a function of cou- 
pling strength /3, evaluated for two different collapse redshifts 
2c = (solid, red) and Zc = 5 (long-dashed, green). 



For small /3 ^ 0.4, (Sc(/3) at Zc > 5 is well described by 
a simple quadratic fitting forinula. 



5^{P) = 1.686(1 -f- al3^) , a ^ 0.556 



(61) 



FIG. 7. Extrapolated linear density contrast at collapse for 
coupled quintessence models with different coupling strength 
/3. For all plots we use a constant a = 0.1. We also depict 5c 
for reference ACDM (dotted, pink) and EdS (double-dashed, 
black) models. 



As opposed to the results found in [52|, no oscillations 
are seen in 5c{zc)- Furthermore, the effect of the cou- 
pling on the extrapolated linear density contrast at col- 
lapse is smaller, though we observe an increase of 5c with 
increasing coupling strength /3, as depicted in Fig|8] for 
two collapse redshifts Zc = and Zc = 5. A coupling 
/3 = corresponds to a ACDM cosmology, hence the ob- 
served 5c is given by 5c ~ 1.686 for z Zc = 5 and by the 
accordingly reduced value for Zc = 0. An increase of /3 
results in an increase of 5c for both redshifts. The reason 



For larger /3 < 1 a fit requires an additional correction 
and reads 



5ciP) = 1.686(1 + a/3^ ~ 5/3^), a = 0.556, h = 0.107 . 

(62) 

It is worth noting that our values of 5c were obtained 
under the assumption that baryonic contributions may be 
neglected, in order to be able to relate the results to the 
simple Einstein de Sitter scenario. Indeed, a numerical 
analysis under inclusion of a baryonic component shows 
a quite significant increase of the critical density contrast 
(5c, leading to values close to those found in [53[. 

Also note that we have limited our analysis to the crit- 
ical density contrast at collapse. Other works 4l| have 
rather focused on the respective quantities at virializa- 
tion. Since there are no fundamental differences, these 
may as easily be evaluated within our formalism. 
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VI. GROWING NEUTRINOS 

Another interesting framework, analogous to coupled 
quintessence, in which a fifth force is present, is the grow- 
ing neutrinos scenario |lll . Il2| . Here, rehc neutrinos ob- 



tain a growing, cosmon dependent mass, implemented by 
a large, negative coupling /3. In this context, dark energy 
domination and the late acceleration of the Universe can 
be naturally explained by relating it to a "trigger event" , 
the recent transition of neutrinos to the nonrelativistic 
regime. 



A. Cosmological model 




1000 



1+z 



As neutrinos have been relativistic particles through 
most of the history of the Universe, Eqs. (|15p . ([T8| and 
([rs]) are appropriately altered to include neutrino pres- 
sure 



iH{l + w^)p^ + (i^{l~iw,)p, (63) 
/3 0(l-3i&,)p, (64) 
(65) 



Py = -3iJ(l + Wy)p, 



As opposed to the models of coupled CDM discussed 
above, the constant /? is now negative and its modulus 
much larger than one. Bounds for the couplings a and j3 
have been discussed in ll|, |67| . For the following analy- 
sis we choose a = 10 and several values for the coupling 
(3 = —52, —112 and —560. Note that the couplings may 
be related to the present neutrino mass via 



TO^(io) 



P 



174to)16cV 



(66) 



where fi^ is the dark energy density fraction today. 

A numerical integration of (p5)) - (p5)) and the appropri- 
ate equations for radiation and CDM leads to the back- 
ground evolution depicted in Fig. |9l While the cosmon 
is on the matter (radiation) attractor at early times, the 
transition of neutrinos to the nonrelativistic regime al- 
most stops the evolution of the cosmon. The dark en- 
ergy density is able to overcome all other components and 
dominates the Universe &i t = t^. As the kinetic contri- 
bution to the cosmon energy density is greatly reduced, 
the latter is dominated by the potential V{(t)), success- 
fully mimicking the behavior of a cosmological constant. 



B. Spherical collapse and growing neutrino 
quintessence 

We have applied our method to the case of growing 
neutrino quintessence. To our knowledge this is the first 
time that spherical collapse is performed on this class of 
models. 



FIG. 9. Energy densities of neutrinos (solid, red), cold dark 
matter (long dashed, green), dark energy (dot-dashed, blue) 
and photons (short dashed, black) are plotted vs redshift. We 
use a sample model with constant P — —52, a = 10 and a 
large average neutrino mass m,^ = 2.11 eV. 



Because of the strong cosmon-mcdiatcd attractive 
force between neutrinos, bound neutrino structures may 
form within these models 68[ . It was shown in [13[ that 
their formation will only start after neutrinos become 
nonrelativistic. A nonlinear treatment of the evolution 
of neutrino densities is thus only required for very late 
times, and one may safely neglect neutrino pressure as 
compared to their density, which substantially simplifies 
the scenario. All calculations of the previous sections are 
thus equally valid for the growing neutrinos scenario; we 
can straightforwardly apply the evolution equations (|52|) 
and (|53p for the nonlinear and linear neutrino density 
contrast. 

In Fig[TO]we plot the evolution of the nonlinear density 
contrast as obtained from numerically solving Eq. (|52p for 
a model with /3 = —52. The linear density contrast, solu- 
tion of Eq. (|53p , is also shown. For comparison, wc have 
also included the linear density contrast resulting from 
the full relativistic equations as given in 13| . The results 
of the linearized spherical collapse and the relativistic 
theory can be seen to agree remarkably well. Compari- 
son with the full hydrodynamic results from Ref. [32| also 
yields a one-to-one agreement, as expected since the lat- 
ter is the basis for the present work. The slight deviation 
around redshift z '^ 1.5 can be accounted for by a short 
recuperation of neutrino pressure at this time. However, 
no significant impact on the extrapolated linear density 
contrast Sc was found. 

In order to illustrate the dependence of the growth of 
the overdensity on the coupling /3, we show the evolution 
of dy{z) and S^^lIz) in Fig. [TT]for the three given cou- 
plings /3 = —52, —112, —560 and a = 10. Given a fixed 
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self-interaction a, a larger (3 corresponds to a smaller 
present neutrino mass; as a consequence, neutrinos be- 
come nonrelativistic at smaller redshifts znr. In the rel- 
ativistic regime, no growth of neutrino perturbations is 
observed in the linear regime. To comply with this, we 
only start the integration of the spherical collapse equa- 
tions once the transition to the nonrelativistic regime is 
observed at redshift 2;nr(/3)- It can be observed in FigfTTJ 
that a higher /3 leads to a strongly enhanced growth of 
the density contrast. On the other hand, because of the 
later transition to the nonrelativistic regime, perturba- 
tions start to grow at much lower redshifts. 
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FIG. 10. Evolution of neutrino nonlinear (solid, red) and 
linear (long-dashed, green) density contrast. For comparison, 
we have also included the relativistic linear density contrast 
including pressure terms (short-dashed, blue). 



The extrapolated linear density at collapse 5c for grow- 
ing neutrino quintessence reflects in all respects the char- 
acteristic features of this model and result in a 5c which 
looks quite different from standard dark energy cosinolo- 
gies. We have plotted the dependence of 5c on the col- 
lapse redshift Zc in Fig jT^] for all three couplings. 

The oscillations seen are the result of the oscillations 
of the neutrino mass caused by the coupling to the scalar 
field: the latter has characteristic oscillations as it ap- 
proaches the minimum of the effective potential in which 
it rolls, given by a combination of the self-interaction po- 
tential U{(j)) and the coupling contribution /3(1 — ?iWy)py. 
Furthermore, due to the strong coupling /?, the average 
value of (5c is found to be substantially higher than 1.686. 
Such an effect can have a strong impact on structure for- 
mation and 5c can then be used within a Press-Schechter 
formalism. 

For the strongly coupled models, corresponding to a 
low present day neutrino mass TOiy(to), the critical density 
at collapse is only available for Zc ^ 0.2, 1 for /3 = —560, 
— 112, respectively. This is again a reflection of the late 
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FIG. 11. Evolution of neutrino nonlinear / linear density 
contrast for a — 10 and fi = —52 (solid, red / long-dashed, 
green), /3 = —112 (short-dashed, blue / dotted, pink) and /3 = 
—560 (dot-dashed, light blue / double-dashed, black). The 
chosen couplings correspond to a present average neutrino 
mass of m^{to) = 2.1 eV, 1 eV and 0.2 eV, respectively. 



transition to the nonrelativistic regime. 

A full nonlinear investigation of single lumps within 
growing neutrino quintessence was performed in [32|. 
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FIG. 12. Extrapolated linear density contrast at collapse Sc 
vs. collapse redshift Zc for growing neutrinos with (3 = —52 
(solid, red), /3 = —112 (long-dashed, green) and (3 = —560 
(short-dashed, blue) . A reference EdS model (double-dashed, 
black) is also shown. 
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VII. EARLY DARK ENERGY 
A. Cosmological model 

A convenient way to parametrize the presence of a non- 
negligible homogenous dark energy component at early 
times (from now on labeled as EDE) was presented in 
[69| . Here, the dark energy density is 



1.688 



PDe(z) == /ODE,0 (1 + z) 



3(l+iDfe(z)) 



with 



PDE,0 = Pcrit.O ^DE,0 = SiJp (1 — ^m.o) 

and the equations of state parametrized by: 

Wo 



Whiz) 



1 



(67) 
(68) 

(69) 



61n(l + z) ' 

where 5 is a constant related to the amount of dark energy 
present at early times 

3iDo 



b=-- 



In 



l-Slr 



In 



i-o„ 



(70) 



Here the subscripts "0" and "e" refer to quantities cal- 
culated today or early times, respectively. Dark energy 
pressure will be given by pde(-z) = Whiz) Pue{z)- If we 
specify the spherical collapse equations for this case, the 
nonlinear evolution of the density contrast follows the 
evolution equations ([5^ and (|53p without the terms re- 
lated to the coupling: 



S„, = -2H5r, 



4 5: 



0711. L 



31 

-2H5„i,L + T^Sm.L A 



^ + -(5„(l + (5„)p™ ,(71) 



(72) 



As before, we assume relativistic components to remain 
homogenous. 



B. Spherical collapse and EDE 

In the following we present our results for two models 
of early dark energy, namely model I and II from [5J]. 
Model I is given by the set of parameters 



i in 



0.332, Wo = -0.93, 17DE,e = 2-10" 



(73) 



whereas model II is parametrized by 

n„i,0 = 0.314 , Wo = -0.99 , f^DE.e 



10"^ . (74) 



Our results for Sc in both EDE models are plotted in 
Fig.IH together with 6c in ACDM. Coherent with the 
results of [56| , we find a suppression of Sc as compared to 
ACDM that is much lower than in the original paper J5J] . 
More precisely, while the latter found Sc{zc = 5) ~ 1.62 
for model I, corresponding to a relative change of ~ 4%, 
we obtain Sc{zc = 5) - 1.685 (~ 5 • 10-2%). 




ACDM 
EDE, model I 
EDE, model II 



FIG. 13. Extrapolated linear density contrast at collapse Sc 
vs. collapse redshift Zc for EDE models 1 (solid, red) and 
II (long-dashed, green), as well as ACDM (double-dashed, 
black) . 



VIII. CONCLUSIONS 

Spherical collapse is a semi analytical method often 
used to estimate the nonlinear evolution of structures 
without reverting to complex numerical methods like N- 
body simulations. 

After reviewing its application to standard cosmolo- 
gies, we have considered the case of coupled dark en- 
ergy cosmologies, in which a fifth force other than gravity 
modifies the collapse. 

We have shown that the inclusion of the fifth force 
within the spherical collapse picture deserves particular 
caution. As spherical collapse is intrinsically based on 
gravitational attraction via the Friedmann equations, it 
does not account for other external forces unless it is 
suitably modified. 

We have presented a detailed comparison between the 
linearized standard spherical collapse picture and the lin- 
ear relativistic equations, whose results are summarized 
in Table HI Applying standard spherical collapse equa- 
tions to coupled dark energy by adding a coupling in the 
conservation equations is insufficient to describe the cou- 
pled dark energy scenario, as the fifth force is still miss- 
ing entirely from the evolution of the density contrast S. 
Results in Table U also show that a standard treatment 
of spherical collapse may lead to problems even in the 
uncoupled case, whenever the scalar field is treated as 
inhomogeneous . 

We have illustrated in detail how a modification of the 
spherical collapse picture which correctly accounts for the 
presence of a fifth force is still possible. We have de- 
rived the set of second order differential equations for the 
density contrast from the fully nonlinear Navier- Stokes 



15 



equations. We have then shown how 5c can be evalu- 
ate directly from these equations and how the spheri- 
cal collapse formalism can be reformulated starting from 
them. Most importantly, we have further checked that 
our results match the numerical resolution of the nonlin- 
ear hydro dynamical equations performed as described in 

Mm- 

We have applied our procedure to coupled quintessence 
scenarios, evaluating the extrapolated linear density at 
collapse for this class of cosmologies and showing how 
it depends on the coupling /?. Furthermore, we have 
for the first time applied the spherical collapse to the 
case of growing neutrino quintessence, where neutrinos 
feel a fifth force interaction that can lead them to clus- 
ter at very large scales. In this case, we demonstrate 
how the extrapolated linear density at collapse shows a 
characteristic oscillating behavior, different from stan- 
dard dark energy models. In future work, this result 
could be used within a Press-Schechter [57| formalism to 
estimate neutrino halo mass distributions. We have fur- 
ther commented on the choice of initial conditions, whose 
choice has to be made with careful attention when dealing 
with the extrapolated linear density at collapse. Finally 
we have used our approach to verify results found in [56| 
on spherical collapse and early dark energy (EDE). 
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Integration yields the following expression for the veloc- 
ity: 



Vm(r, t) 



c 



47rr2(l + 5,„) 



Therefore 



lim ( —^ ] 
r-yO \ r J 



e 



r™ y 47rr3(l + 6,n) 



(78) 



(79) 



Expanding S{r',t) in ([77)) around r' = enables us to 
rewrite ^ as 

ar,t) = 47r / r'2 {6„,{0,t) + 0{r')) dr' 
Jo 
4-7r 
= — r3<5„,(0,t) + O(r4) 



(80) 




rA^6U^r3(l + <5™) ) ~ 3^-^'"lx=o 



We obtain 



V(v„V)v„4^o = - (V-v„0'|,^o (81) 
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APPENDIX 

In order to derive the identity ()51l it is convenient to 
express V(vm V)vm in spherical coordinates. Symmetry 
implies v„i = Vm e,. and therefore 



V(v„ V)v„i = 



(V • v„ 



f V„ V • V„, + V„j d^Ym 



(75) 

which we may evaluate at x = making use of (|47p and 
0: 



2V • v„i lim 



r— ^0 T 



V(v™ V)v™|^^o = (V • v„)' 

■/■ — ?u /■ 

(76) 

To evaluate the limit, we consider a scalar function de- 
fined as 



i{r,t):=ATT / r'^5,n{r',t)dr' 



(77) 



We differentiate it in time and insert (|^^ in spherical co- 
ordinates, using V- [(1 -I- 5m)^rn] = J^dr ^(1 "f <5m)v„] . 
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